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1.  Introduction 


The  prime  function  of  the  respiratory  system  is  to  exchange  oxygen  from  the  air  into 
the  venous  blood  and  carbon  dioxide  from  venous  blood  to  air.  Gases  are  transported 
between  the  environment  to  the  alveoli  by  convection.  The  pulmonary  airways  bifurcate  in 
a  treelike  structure  as  shown  in  Fig.  1  [16].  Hence  the  geometry  of  the  airways,  such  as 
diameters  and  lengths,  airway  wall  thickness,  and  branching  patterns  are  essential  physical 
properties  of  the  system  that  significantly  influence  the  transport  of  respiratory  gases. 


.  STRUCTURE  AND  FUNCTION 


Figure  1 .  Cast  of  the  airways  of  a  human  lung.  The  alveoli  have  been  pruned  away, 
but  the  conducting  airways  from  the  trachea  to  the  terminal  bronchioles  can  be  seen  [16]. 
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The  forced  oscillation  technique  is  a  widely  used  approach  to  extract  insight  into  how 
the  structure  of  the  lung  (airway  geometry  and  mechanical  properties)  influences  it’s 
function.  There  are  two  common  methods  for  applying  forced  oscillations  to  measure 
mechanical  impedance.  One  is  to  impose  forced  oscillations  at  the  airway  opening  and 

measure  the  complex  ratio  of  airway  opening  pressure  Pao  to  flow  V  «0 ,  which  yields  input 

impedance  [6,  25],  The  other  is  to  apply  forced  oscillations  around  the  chest 

wall  Pew  and  measure  transfer  impedance  defined  as  ZU=PCW/ V  (Fig.  2).  In  either  case, 
the  impedance  spectrum  is  analyzed  to  understand  the  mechanical  properties  of  the 
respiratory  system.  Conceptually,  the  respiratory  system  is  comprised  of  airways  and 
tissues  separated  by  a  volume  of  compressible  alveolar  gas.  The  mechanical  properties 
and  morphometric  structure  of  the  airways  and  tissues  produce  effective  resistances, 
inertances  and  compliances,  which  all  effect  the  overall  impedance  spectrum. 
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(Oi  > 


Transfer  Impedance  System 


(k  )  |  Input  Impedance  Systems 


Standard  Method:  Zst 


-  Experimental  set-up.  The  conscious  mouse,  placed  in  the 
modified  Battelle  tube  in  the  body  box,  respires  from  the  front  cham¬ 
ber.  To  prevent  rcbreathing,  the  from  chamber  is  flushed  with  a  bias 
flow.  R:  wire  mesh  screen  resistance:  Pbs:  pressure  in  the  body  box 
(body  surface) ;  Pao:  pressure  in  the  from  chamber  (airway  opening). 


Figure  2.  (a)  and  (b)  Schematic  of  measurement  systems  for  input  and  transfer 
impedances  in  humans,  respectively  [6],  (c)  Schematic  of  measurement  system 
for  transfer  impedance  in  rats  [25]. 
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A  theoretical  study  of  a  dog  lung  model  predicted  that  with  forced  oscillations 
imposed  at  the  airway  opening,  the  alveolar  pressure  Pa|V  could  be  higher  than  the  applied 
pressure  Paoat  low  frequencies  (below  8  Hz)  [20],  The  question  remained  whether  similar 
pressure  amplification  would  occur  in  the  human  lung.  Furthermore,  would  there  be 
pressure  amplification  within  the  airways  when  the  oscillations  are  applied  to  the  body 
surface?  Should  there  be  any  pressure  amplification,  what  would  be  the  implication  for 
patients  with  bronchoconstriction? 

The  main  goals  of  this  study  were  to  investigate  the  lung  mechanics,  including  input 
impedance  (Zin),  transfer  impedance  (Z*),  as  well  as  pressure,  strain,  and  wave 
propagation  distributions,  in  order  to  answer  these  questions. 

The  computational  predictions  were  based  on  the  symmetrical  and  asymmetrical 
models  of  the  human  lung.  Implementation  of  the  symmetrical  model  is  relatively 
straightforward  but  the  resulting  predictions  may  not  be  entirely  appropriate  since  the  real 
airways  do  not  branch  symmetrically.  More  accurate  predictions  may  be  obtained  with  an 
asymmetrically  branching  airway  network,  such  as  the  Horsefield  model.  Implementation 
of  an  asymmetrical  model  to  predict  input  impedance  is  also  rather  straightforward  but 
one  to  predict  transfer  impedance  is  very  difficult.  To  solve  the  problems  arising  from 
modeling  an  asymmetrical  tree  structure,  a  superposition  model  was  proposed  based  on 
the  assumption  that  the  system  behaves  linearly  when  small  amplitude  forced  oscillations 
are  applied.  The  tree  structure  of  the  airway  was  first  horizontally  divided  into  an 
asymmetrical  zone  and  a  symmetrical  zone.  The  asymmetrical  zone  was  further  pruned 
into  three  subtrees.  A  uniquely  indexed  binary  tree  representation  was  developed  from 
Horsefield’s  model  using  recursive  algorithms  applied  to  the  asymmetrical  zone  to 
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compute  the  flow  in  all  airway  segments  resulting  from  pressure  sources  at  all  alveoli, 
each  considered  separately.  The  sum  of  these  flows  was  then  calculated  and  the  transfer 
impedance  was  computed  straightforwardly  from  its  definition. 

A  rat  lung  model  was  also  implemented  to  aid  in  the  design  of  experiments  to 
determine  sound  exposure  limits  in  humans. 
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2.  Background 


2.1.  Symmetrical  airway  morphometry:  Weibel’s  model  of  human 
lung  and  Yeh’s  model  of  rat  lung 

The  treelike  structure  (Fig.  1)  of  the  pulmonary  airways  can  be  analyzed  in 
mathematical  terms  and  related  to  its  function  [16],  Two  commonly  used  human 
airway  models  are  the  Weibel’s  symmetrically  branching  tree  structure  [11]  and  the 
Horsefield’s  asymmetrically  branching  tree  structure  [12].  In  the  human  lung  model 
proposed  by  Weibel,  the  airways  branch  24  times  from  the  trachea  (order  1)  to  the 
alveolar  ducts  (order  24).  Each  branch  divides  into  two  identical  daughter  branches. 
The  lengths  and  diameters  of  the  airways  were  determined  by  standard  morphological 
techniques  on  five  normal  human  lungs  inflated  to  total  lung  capacity  (TLC)  [11].  All 
airways  of  a  given  order  are  assumed  to  have  the  same  diameter  and  length.  Thus,  all 
pathways  from  the  trachea  to  the  alveoli  are  identical.  The  mean  diameter  and  length 
of  each  airway  order  are  listed  in  Table  1. 
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Table  1.  Weibel’s  symmetrical  model  A  of  the  human  lung 


Order  n 


Number  per  generation  TV 


64 


128 

256 


512 


1024 


2048 


4096 


8192 


16384 

32768 


65536 


131072 


262144 


524288 


1048576 


2097152 


4194304 


8388608 


Diameter  d(n)  cm 

Length  l(n)  cm 

1.8 

12 

1.221 

4.78 

0.828 

1.91 

0.562 

1.27 

0.445 

1.07 

0.351 

0.9 

0.281 

0.76 

0.227 

0.76 

0.186 

0.64 

0.154 

0.54 

0.13 

0.46 

0.11 

0.39 

0.095 

0.33 

0.083 

0.27 

0.073 

0.23 

0.065 

0.2 

0.059 

0.16 

0.054 

0.14 

0.05 

0.12 

0.047 

0.1 

0.045 

0.08 

0.043 

0.07 

0.042 

0.06 

0.042 

0.05 

Unlike  Weibel’s  model,  the  model  of  the  rat  lung  proposed  by  Yeh  et  al  [16]  is  a 
model  without  connectivity.  Based  on  very  complete  anatomical  data  measured  by 
Raabe  et  al.,  the  number  of  terminal  bronchioles  (Nth)  was  measured  and  the  median 
number  of  generation  (order)  m  down  all  the  pathways  to  terminal  bronchioles  was 

•  m~l 

calculated.  The  average  branching  ratio  x  was  obtained  from  the  equation  No,  =  x 
where  x  may  not  be  an  integer,  which  implies  that  the  connectivity  is  undefined.  This 
may  also  be  called  an  “average-path”  lung  model.  The  average-path  model  of  the 
airways  in  the  rat  lung  is  shown  in  Table  2.  To  be  consistent  with  the  numbering  used 
in  our  Weibel  model,  order  24  in  Table  2  represents  the  trachea,  and  generation  1 
represents  the  terminal  airway. 
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Table  2.  Yeh’s  symmetrical  model  of  rat  lung 


order  n 

Number  per  generation  N 

Diameter  d(n)  cm 

Length  l(n)  cm 

24 

1 

0.3400 

2.6800 

23 

2 

0.2900 

0.7150 

22 

3 

0.2630 

0.4000 

21 

5 

0.2030 

0.1760 

20 

8 

0.1630 

0.2080 

19  j 

14 

0.1340 

0.1170 

18 

23 

0.1230 

0.1140 

17 

38 

0.1120 

0.1300 

16 

65 

0.0950 

0.0990 

15 

109 

0.0870 

0.0910 

14 

184 

0.0780 

0.0960 

13 

309 

0.0700 

0.0730 

12  ^ 

521 

0.0580 

0.0750 

11 

877 

0.0490 

0.0600 

10 

1477 

0.0360 

0.0550 

9 

2487 

0.0200 

0.0350 

8 

4974 

0.0170 

0.0290 

7 

9948 

0.0160 

0.0250 

6 

19896 

0.0150 

0.0220 

5 

39792 

0.0140 

0.0200 

4 

79584 

0.0140 

0.0190 

3 

159168 

0.0140 

0.0180 

2 

318336 

0.0140 

0.0170 

1 

636672 

0.0140 

0.0170 

2.2.  Asymmetrical  airway  morphometry:  Horsefield’s  model  of  human 
lung 

In  the  Horsefield  model  the  airways  are  classified  by  the  airway  order  n  starting 
from  the  terminal  airway  (order  1)  to  the  trachea  (order  36)  [12].  In  this  model,  an 
airway  of  order  n  branches  into  two  daughters  of  order  n-1  and  order  n-l-A(n),  where 
A(n)  is  the  recursion  index  that  characterizes  the  degree  of  asymmetry  at  each  airway 
order.  As  one  can  see,  connectivity  in  this  model  is  defined.  The  model  represents  the 
average  asymmetry  found  in  the  lung  studied  by  Horsefield  et  al.  The  diameter  d(n), 
length  l(n),  as  well  as  a  recursion  index  A(n)  of  each  airway  order  are  specified  (Table 

3)- 

Note  that  we  have  modified  the  diameters  of  the  last  eight  orders  since 
physiologically,  it  is  unlikely  to  have  increased  diameters  for  lower  orders.  The  new 
diameters  were  assumed  to  follow  the  diameter-generation  relationship  given  by  Weibel 
[11]: 


d(n)  =  do  •  e-(0'293-000624-n)-n 

where  do  is  diameter  of  9th  order  in  Horsefield’s  model  which  is  the  starting  point  of  the 
fitting  curve,  and  n  is  the  generation  number  in  Weibel’s  model. 

The  comparison  of  modified  diameters  and  the  original  ones  are  shown  in  Fig.  3. 
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Table  3.  Horsefield’s  asymmetrical  model  of  human  lung 


Order  n 

Number  per  order# 

Delta  A(n) 

Diameter  d(n)  cm 

Length  l(n)  cm 

36 

1 

1 

1.6000 

10.000 

35 

1 

2 

1.2000 

5.000 

34 

2 

3 

1.1000 

2.200 

33 

2 

3 

0.8000 

1.100 

32 

3 

3 

0.7300 

1.050 

31 

3 

3 

0.5900 

1.130 

30 

5 

3 

0.5900 

1.130 

29 

7 

3 

0.5400 

0.970 

28 

10 

3 

0.4300 

1.080 

27 

13 

3 

0.3500 

0.950 

26 

18 

3 

0.3500 

0.860 

25 

25 

3 

0.3100 

0.990 

24 

35 

3 

0.800 

23 

48 

3 

0.2800 

0.920 

22 

66 

3 

0.2700 

0.820 

21 

91 

3 

0.810 

20 

126 

3 

0.2400 

0.770 

19 

174 

3 

0.2180 

0.640 

18 

240 

3 

0.2000 

0.630 

17 

331 

3 

0.1800 

0.517 

16 

457 

3 

0.1600 

0.480 

15 

631 

2 

0.1400 

0.420 

14 

871 

2 

0.1100 

0.360 

13 

1202 

1 

0.310 

12 

2290 

0 

0.0760 

0.250 

11 

6653 

0 

0.0630 

0.110 

10 

13306 

0 

0.0530 

0.131 

9 

26612 

0 

0.0480 

0.105 

8 

53224 

0 

0.0430 

0.075 

7 

106448 

0 

0.0323 

0.059 

6 

212896 

0 

0.0245 

0.048 

5 

425792 

0 

0.0189 

0.048 

4 

851584 

0 

0.0147 

0.048 

3 

1703168 

0 

0.0116 

0.048 

2 

3406336 

0 

0.0093 

0.048 

1 

6812672 

0 

0.0075 

0.048 
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3.  Methods 


3.1  Transfer  function  for  single  airway  segment 

3.1.1  Series  impedance  and  shunt  admittance  of  a  nonrigid  tube 

The  acoustical  properties  of  an  airway  segment  were  modeled  in  terms  of  the  series 
impedance  composed  of  gas  resistance  (Rg)  and  inertance  (Ig),  and  shunt  admittance  which 
is  due  to  gas  compression  (Cg),  thermal  losses  (Gg),  and  wall  properties  (Zw)  (Fig.  4) 
[3,15]. 


Figure  4.  Transmission  line  representation  of  small  segment  of  nonrigid  tube 


The  series  impedance  Z&s  and  shunt  admittance  (Yg>p+Yw)  are  defined  as 


(i) 

,  +  Yw  =  Gg  +jcoCg+Yw 

(2) 
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where  Yw  is  wall  admittance  (Yw  =  1/  Zw).  The  coupling  of  Zg,s  and  shunt  admittance 
determines  the  characteristic  impedance  Z0  and  propagation  coefficient  y  given  as 


Zo  = 


g.s 

+  Yw 


y  =  -y/Zg.s  •  (Yg.  p  +  Yw)  =  a  +  jp 


(3) 

(4) 


where  a  is  the  attenuation  coefficient  determining  the  amount  of  exponential  attenuation  of 
the  wave  amplitude  along  the  airway  segment;  and  (3  is  the  measure  of  phase  change  of  the 
wave.  The  wave  propagation  velocity  c  and  /?  are  related  by 


© 

c  =  — 

P 

where  ©  =  2 %f,  /is  frequency  in  Hz. 


(5) 


3.1.2  Transmission  line  theory  and  two-port  network  representation 

Each  segment  of  the  airway,  which  is  the  basic  computational  unit  of  the  proposed 
model,  can  be  modeled  as  a  two-port  (four-terminal)  network  as  illustrated  in  Fig.  5. 
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Pa,  Qa 


Figure  5.  Four-terminal  network  representation  of  a  segment  of  airway 
The  pressure  drop  and  flow  discharge  of  the  segment  based  on  transmission  line 
theory  are  given  by  [10] 

Pu  =  Pj  cosh(y/)  +  QdZo  sinh(y/) 

Qu  =  Pd  sink(Y0.+  Qd  cosh(y/) 

JL.Q 

where  Pu  and  Pa  are  the  upstream  and  downstream  pressure  (in  cmF^O),  respectively.  Qu 
and  Qd  are  the  upstream  and  downstream  flow  (in  L/s),  respectively.  And  /  is  the  length  of 
the  airway  segment  (in  cm).  The  characteristic  impedance  Zo  and  propagation  coefficient  y 
are  obtained  from  Eq.  (3)  and  (4). 

The  transmission  matrix  (T-parameter)  for  the  network  in  Fig.  5  can  be  derived  from 
Eq.  (6)  and  Eq.  (7)  as 


(6) 

(7) 


Pu,  Qu 


~Pu~ 

~A 

B~ 

~Pd~ 

_Q«_ 

C 

D_ 

_Q*_ 

(8) 
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where 


A  =  cosh(?/) 

(9) 

5  =  Z0sinh(^) 

(10) 

c  _  sinh(^) 

01) 

Zo 

D  =  cosh(  yi  ) 

(12) 

In  the  implementation  of  the  model,  these  airway  segments  (equivalent  networks)  are 
connected  either  in  series  or  in  parallel.  Connecting  segments  in  series  is  conveniently 
done  by  multiplying  T-matrices  of  all  the  network  blocks.  However,  connecting  segments 
in  parallel  is  more  conveniently  accomplished  by  adding  their  admittances,  Y-parameters, 
defined  by 


Qu 

Qa 


cosh(yl) 

1 

yn  yi2 

"Pu" 

Zo  sinh(yl) 

Zo  sinh(yl) 

y2i  y22_ 

Pd 

1 

cosh(yl) 

-  J 

Zo  sinh(yl) 

Zosinh(yl) 

Qu 

Qd 


(13) 


3.1.3.  Model  of  airway  wall  properties  and  tissue  properties 

The  measurement  of  wave  propagation  velocity  in  dog  tracheae  made  by  Guelke  and 
Bunn  [21]  indicated  that  it  was  frequency  dependent  which  implicates  that  the  airway  is 
not  rigid.  A  single  compartment  model  including  resistive,  inertive,  and  elastic  elements 
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was  used  to  explain  this  frequency  dependence.  Later,  Suki  et  al.  [15]  found  a  more 
complex  frequency  dependence  of  propagation  velocity  in  excised  calf  tracheae  which  led 
to  a  two-compartment  representation  of  wall  properties  [15].  As  suggested  by  their  study, 
a  two-compartment  model  is  used  to  model  soft  tissue  and  cartilage  in  the  airway  walls 
with  diameters  greater  than  2  mm  (Fig.  6).  The  walls  of  airways  whose  diameter  is  less 
than  2  mm  are  modeled  with  a  single  soft  tissue  compartment  due  to  the  absence  of 
cartilage. 


Soft  Tissue 


Cartilage 


Figure  6.  Two  compartment  model  of  airway  wall  properties 
The  volumetric  properties  of  the  airway  walls  are  computed  from  equations  derived  by 
Suki  et  al.  [15]  and  given  below.  The  resistance  is  given  by 


Rw(n)  = 

7td3(n)l(n) 


(14) 


where  d(n),  l(n),  and  h(n)  are  airway  segment  diameter,  length  and  wall  thickness, 
respectively;  and  v  is  viscosity  of  either  cartilage  or  soft  tissue.  The  inertance  is  given  by 


_  h(n)p 

7id(n)l(n) 


(15) 
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where  p  is  density  of  the  tissue  (either  cartilage  or  soft  tissue).  Finally,  the  capacitance  is 
given  by 


Cw(n)  — 


;td3(n)l(n) 

4h(n)Y 


(16) 


where  Y  is  Young’s  modulus  of  either  cartilage  or  soft  tissue. 

The  densities  of  the  cartilage  and  soft  tissue  (pc  and  ps)  were  fixed  to  1.14  and  1.06 
g/ml,  respectively  [1],  The  Young’s  modulus  for  cartilage  (Yc)  was  fixed  to  the  mean 
value  (4000  cmH20)  reported  in  human  tracheal  cartilage  stripes  [19].  The  soft  tissue 
viscosity  vs  and  the  Young’s  modulus  Ys  were  fixed  to  the  mean  values  of  1.04  cmH20 
and  593  cmH20  respectively,  which  were  reported  by  Habib  et  al.  [1].  Based  on  the 
empirical  findings  of  Suki  et  al.  [15],  it  was  assumed  that  cartilage  and  soft  tissue  have 
comparable  viscosity-to-elasticity  ratios;  i.e.,  Vs/Ys  =  Vc/Yc .  Therefore,  the  cartilage 
viscosity  vcwas  fixed  to  a  value  of  (7.015  cmH20)  calculated  from  this  formula. 

The  cartilage  content  c(n)  is  a  function  of  d(n)  as  proposed  by  Gunst  and  Stropp  [15]: 


c(n)  =  -0.09d2  (n)  +  0.25d(n)  -  0.024  (17) 


Since  the  soft  tissue  content  is  simply  1  -c(n),  the  effective  wall  volumetric 
impedance  Zw,efr  of  an  airway  segment  of  order  n  is 


Zw,  cfr(n) 


Zw,  c(n)  •  Zw,  s(n) 

c(n)  •  Zw,  s(n)  +  [1  -  c(n)] "  Zw.  c(n) 


(18) 
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where  Zw,c  and  Zw>s  represent  cartilage  and  soft  tissue  compartment  impedances, 
respectively. 

The  wall  thickness  is  required  to  compute  the  volumetric  properties  of  the  airway 
wall  (Eqns.  14  -  Eq.  16).  Habib  et  al.  [1]  used  a  polynomial  expression  approximating  the 
data  of  the  wall  area  internal  to  the  smooth  muscle  layer  WAsm(n)  as  a  function  of  internal 
airway  perimeter  7td(n)  for  healthy  human  lungs  from  Wiggs  et  al.  [15].  By  assuming  a 
circular  geometry  for  both  the  luminal  and  external  wall  surfaces,  the  wall  thickness 
internal  to  the  outer  smooth  muscle  layer  h^n)  was  computed  as  a  function  of  internal 
diameter  d(n)  as  [1] 

hsm(n)  =  -Jd2(n)  /  4  +  WAsm(n)  /  k  -  d(n)  /  2  09) 

Furthermore,  in  order  to  obtain  the  thickness  of  the  entire  airway  wall  (not  only 
smooth  muscle  layer),  Habib  et  al.  applied  system  identification  methods  to  estimate  the 
total  wall  thickness  as  a  function  of  the  diameter  of  a  given  generation.  In  the  current 
study,  the  airway  thickness  data  of  Habib  et  al.  are  used  for  our  Horsefield-based, 
asymmetrical  human  lung  model.  The  thickness  of  the  airway  walls  in  our  Weibel-based, 
symmetrical  human  model  and  our  Yeh-based  rat  lung  model  were  predicted  from  a 
polynomial  fit  of  thickness  as  a  function  of  airway  order  published  by  Habib  et  al.  [1] 
(Fig-  7). 
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Figure.  7.  Polynomial  fit  of  wall  thickness  data  from  Habib  et  al. 


The  wall  thickness  and  cartilage  content  data  are  listed  in  Table  4,  Table  5,  and  Table 
6,  for  the  Horsefield’s  human  lung  model,  the  Weibel’s  human  lung  model,  and  the  Yeh’s 
rat  lung  model,  respectively. 
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Table  4.  Airway  wall  parameters  based  on  diameters  of  the  Horsefield’s  model 


Order  n 

Wall  thickness  h(n)  cm 

Cartilage  Content  c(n) 

36 

0.6554 

0.6700 

35 

0.3054  J 

0.5000 

34 

0.2372 

0.5000 

33 

0.0929 

0.3300 

32 

0.0720 

0.2500 

31 

0.0429 

0.2000 

30 

0.0429 

0.0926 

29 

0.0361 

0.0851 

28 

0.0262 

0.0671 

27 

0.0222 

0.0526 

26 

0.0222 

0.0526 

25 

0.0208 

0.0450 

24 

0.0201 

0.0410 

23 

0.0197 

0.0390 

22 

0.0195 

0.0370 

21 

0.0188 

0.0329 

20 

0.0185 

0.0309 

19 

0.0176 

0.0262 

18 

0.0169 

0.0224 

17 

0.0160 

0 

16 

0.0150 

0 

15 

0.0139 

0 

14 

0.0118 

0 

13 

0.0106 

0 

12 

0.0088 

0 

11 

0.0074 

0 

10 

0.0063 

0 

9 

0.0056 

0 

8 

0.0051 

0 

7 

0.0038 

0 

6 

0.0029 

0 

5 

0.0022 

0 

4 

0.0017 

0 

3 

0.0014 

0 

2 

0.0011 

0 

1 

0.0009 

0 

! 
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Table  5.  Wall  parameters  for  Weibel’s  model. 


Order  n 

h(n) 

c(n) 

24 

0.9321 

0.1344 

23 

0.3059 

0.1471 

22 

0.1094 

0.1213 

21 

0.0460 

0.0881 

20 

0.0303 

0.0694 

19 

0.0215 

0.0527 

18 

0.0167 

0.0391 

17 

0.0137 

0.0281 

16 

0.0118 

0.0194 

15 

0.0104 

0.0124 

14 

0.0095 

0.0070 

13 

0.0087 

0.0024 

12 

0.0082 

0 

11 

0.0077 

0 

10 

0.0074 

0 

9 

0.0071 

0 

8 

0.0069 

0 

7 

0.0067 

0 

6 

0.0066 

0 

5 

0.0064 

0 

4 

0.0064 

0 

3 

0.0063 

0 

2 

0.0063 

0 

1 

0.0063 

0 
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Table  6.  Wall  parameters  for  Yeh’s  model. 
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The  most  peripheral  airways  were  assumed  to  be  uniformly  terminated  by  a  single 
alveolar  sac  (Fig.  8).  The  terminal  parenchymal  and  chest  wall  tissues  are  modeled  as 
series  resistance-inertance-capacitance  tissue  impedance  in  parallel  with  a  capacitance  due 
to  alveolar  gas  compression  [1,5]  (Fig.  9). 


Figure  8.  The  most  peripheral  airway  terminated  by  a  single  alveolar  sac. 


last  airway  segment 


Rti 

A/W 


Cti 


31 


Figure  9.  Model  of  tissue  properties  and  alveolar  gas  compression. 
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The  parameters  (Ru,  Iti,  and  Cti)  of  tissue  properties  of  both  human  lung  models  were 
taken  from  Peslin  et  al.  [8],  The  rat  lung  tissue  properties  were  assigned  the  values 
published  by  Jackson  et  al.  [23],  The  alveolar  gas  compression  compliance  Cg  was 
calculated  with  the  assumption  of  isothermal  conditions  from  assigned  lung  volume  VL 
(usually  takes  the  value  of  functional  residual  capacity  FRC)  and  mean  pressure.  The  gas 

volume  in  a  single  alveolar  sac  is  — ,  where  N  is  the  total  number  of  terminal  airways. 

The  partial  pressure  of  the  compressible  100%  humidified  gas  in  the  alveoli  was  assumed 
to  be  Patm-Prno,  where  Patm  is  atmospheric  pressure  (1033  cmH20)  and  P^ois  the  partial 
pressure  of  water  vapor  at  100%  saturation  (64  CIT1H2O)  [4],  Thus,  Cg  is 

r-  - _  (20) 

N-( Patm-PHio) 


Using  dimensions  (cross-sectional  area  versus  distance  from  the  airway  opening)  of  the 
oral  and  oral  pharyngeal  cavities  (Dan  Olsen,  unpublished  measurements)  we  approximated 
the  geometry  of  the  upper  airways  as  6-segments  whose  lengths  and  diameters  are  given  in 
Table  9,  where  the  wall  thickness  was  computed  from  Eq.  19.  Habib  provided  evidence  that 
the  upper  airways  in  the  human  lung  behave  as  a  two  compartment  model  whose  lumped 
properties  are  given  in  Table  8  [26].  We  know  of  no  data  by  which  we  could  assign 
separate  values  to  the  6  separate  segments  of  our  model.  Therefore  we  simply  assumed  that 
those  airway  wall  properties  suggested  by  Habib  were  distributed  uniformly  along  the 
upper  airways.  Order  1  is  the  segment  next  to  trachea  and  order  6  is  the  uppermost  one. 
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Table  7.  Structural  parameters  of  the  upper  airways  of  the  human  airway  model 


Order  n 

Diameter  d(n)  (cm) 

Length  l(n)  (cm) 

Thickness  h(n)  (cm) 

6 

4.4000 

2.7500 

1.8000 

5 

4.9000 

2.7500 

1.8000 

4 

4.4000 

2.7500 

1.8000 

3 

3.5000 

4.4500 

1.8000 

2 

3.0000 

5.0000 

3.4000 

1 

2.2000 

5.5000 

3.4000 

Table  8.  Lumped  upper  airway  wall  parameters  from  two-compartment  model 


Subject  No. 

Rch 

Ich 

Cch 

Rsm 

k 

CmH^O-s 

CmFLO-s 

L/CmH2 

CmELO-s/L 

1 

14.5 

0.0173 

0.0009 

1.7 

1.18 

2 

75.8 

0.0195 

0.0020 

10.9 

1.7 

3 

9.7 

0.0197 

0.0009 

33.1 

2.23 

4 

33.0 

0.0343 

0.0005 

21.0 

2.75 

5 

56.9 

0.0104 

0.0014 

9.5 

3.53 

Mean 

38.0 

0.0202 

0.0011 

15.2 

2.28 

±SD 

±25.2 

0.0078 

0.0005 

10.8 

0.82 

The  excised  lung  used  in  Weibel’s  study  was  fixed  to  2400-ml  [18].  It  was  assumed 
that  the  lung  expansion  was  isotropic  and  the  diameters  were  scaled  linearly  (i.e.,  70%) 
[24].  Airways  were  assumed  to  be  filled  with  alveolar  gas,  i.e.  100%  humidified  room  air. 
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3.2  Input  impedance  and  transfer  impedance  for  symmetrically  branching  tree 
structure 

3.2.1  Structure  simplification 

Consider  the  symmetrically  branched  tree  structure  (Fig.  10).  The  basic 
computational  unit  of  the  proposed  model  is  a  segment  of  airway,  which  is  modeled  as  a 
two-port  (four-terminal)  network  as  discussed  in  section  3.1.  Hence  the  entire  tree- 
structured  airway  can  be  modeled  as  a  network  illustrated  in  Fig.  11. 


t  t  _ : t _ i  t 

P  (Input  from  alveoli) 


Figure  10.  Symmetrical  airway  tree  structure  with  input  imposed  at  airway  opening  (single 
flow  source)  OR  input  imposed  at  alveoli  (multiple  pressure  sources) 


The  pressure  drops  across  all  airway  segments  of  a  given  generation  are  identical, 
since  the  impedances  of  all  segments  are  identical  (Fig.  1 1).  Hence,  the  pressures  at  points 
a]  and  ci2  are  the  same,  and  pressures  at  points  bj  and  b^  are  the  same  as  well.  For  the  same 
reason,  the  pressures  at  points  C1-C4  are  identical  and  the  pressures  at  points  di-d4  are  also 
identical.  By  connecting  these  equal-pressure  points,  the  simplified  network  structure  is 
obtained  as  shown  in  Fig.  12  [13]. 


Upper 

Airways 


Figure  12.  Simplified  network  structure.  The  dichotomous  branching  pattern  turned 
into  simpler  parallel  connection  of  two-port  network  blocks  for  the 
corresponding  airway  segments. 
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The  structure  in  Fig.  12  can  be  further  simplified  by  substituting  an  equivalent  two- 
port  network  for  each  airway  generation,  namely,  collapsing  all  the  network  blocks  into  a 
single  equivalent  two-port  network  block  with  the  algorithm  described  in  section  3.1.2.  The 
further  simplified  series  two-port  network  structure  is  shown  in  Fig.  13. 

Upper 
Airways 

Jj p* 

nr 

U~«.| 


Trachea 


Figure  13.  Series  two-port  network  equivalent  model  of  symmetrical  airway  tree  structure. 

3.2.2  Algorithm  for  computing  the  equivalent  T-parameters  of  parallel  network 
connection 

The  Y  and  T  parameters  of  one  two-port  network  block  of  a  given  generation  can  be 
obtained  as  described  in  section  3.1.2.  The  equivalent  Y  parameter  matrix  Yp  of  m  blocks 
connected  in  parallel  as  shown  in  Fig.  14  is  [17] 

y„=| y,  (2i) 

1=1 

where  Yi  is  the  Y  matrix  of  /th  network  block,  z=l,  2, . . .  m. 
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Figure  14.  m  two-port  network  are  connected  in  parallel 
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For  m  identical  network  blocks,  where  YpY  for  i=l,  2,  ...  m,  the  parameters  of  the 
equivalent  Yp-matrix  are  simply  the  product  of  m  and  the  corresponding  parameters  in  the  Y- 
matrix: 


YP 


Y  = 

m-yn 

m-yn 

m-yi  i 

m-yn 

(22) 


The  T  parameters  expressed  in  terms  of  the  Y-parameters  are  (see  details  in  Appendix  A) 


'tn 

1 1 2 

y22 

y2i 

, - 

1  >> 

1 

_t21 

t22_ 

<n 

<  |  >> 

I 

_ i 

_  yn 

y2i 

(23) 


where  Ay  is  the  determinant  of  the  Y  matrix.  The  determinant  of  the  Yp-matrix  expressed 


in  terms  of  the  Y-matrix  is 


Ayp  =  m2  •  Ay 


(24) 


Therefore,  Tp  can  be  derived  in  terms  of  the  T-parameters  from  Eq.  (22)  and  Eq.  (23)  as 
follows: 


>> 

e 

i _ 

i 

- 

tl2 

m  •  y2i 

m  •  y2i 
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m 

m2  •  Ay 

m  •  yn 

m  •  1 2i 

t22 

m  •  y2i 

m  •  y2i 

(25) 
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Thus,  even  for  parallel  connections  of  the  network,  we  do  not  have  to  compute  the  Y- 
parameters.  Only  the  T-parameters  are  needed  to  compute  impedance  of  individual  airway 
segments  and  equivalent  networks  of  parallel  connections  of  these  segments. 

In  WeibeFs  model  of  the  human  lung,  m  of  a  given  generation  is  m=  2”,  where  n  is  the 
generation  number.  In  the  Yeh’s  model  of  the  rat  lung,  m  is  calculated  based  on  the  average 
branching  ratio  (which  is  not  2,  and  sometimes,  is  not  even  an  integer)  of  the  given 
generation  as  listed  in  Table  2.  Here,  it  can  be  noticed  that  m  is  not  two  to  the  power  of  the 
generation  number  n. 

3.2.3  Algorithm  for  computing  the  equivalent  T  parameters  of  series  network 
connections 

After  obtaining  the  Tp’s  for  all  generations,  the  next  task  was  to  compute  the 
equivalent  T  parameter  matrix,  Ts>  of  the  series  network  connection  (Fig.  15)  in  order  to  get 
a  single  equivalent  network. 


Figure  15.  Series  connection  of  A:  two-port  network  blocks 
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The  equivalent  T-parameter  matrix  Ts  of  k  network  blocks  connected  in  series  as 


shown  in  Fig.  15  is  [17] 


Ts^YlTp,  (26) 


This  can  be  realized  using  the  following  recursive  algorithm: 
Initialize  Ts  as  a  2x2  unity  matrix, 
for  i=l  to  k 

Ts=  Ts*Tpi 

end 

The  resulting  single  two-port  network  block  is  shown  in  Fig.  16. 


Airway 


Figure  16.  Single  two-port  network  equivalent  of  airway  tree  structure 

3.2.4  Input  and  transfer  impedances 

The  model  for  the  entire  respiratory  system  (Fig.  17)  is  obtained  by  combining  the 
airway  and  tissue  properties. 
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Airway 


Tissue 


Figure  17.  Integrated  respiratoiy  system  model 


For  Zin  computations,  the  impedance  of  the  chest  wall  and  tissues,  Rti-Iti-Cti  in 
series,  which  is  placed  in  parallel  with  alveolar  gas  compression  capacitance,  can  be 
considered  as  the  load  impedance  of  the  airway  two-port  network  (Fig.  17).  Thus  Zin 
looking  from  the  trachea  is  computed  as 


Zin  — 


tsil  *  Zload  +  ts\2 
tszi  •  Zload  +  tsn 


(27) 


where  ts„,  tSn,  ts21,  and  tSM  are  the  four  parameters  in  Ts  matrix  ;  and  Z|0ad  is 


Zload  — 


Zti  •  Zg 
Zti  +  Zg 


(28) 
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where 


Airway  Tissue 


Figure  18.  Integrated  network  model  for  transfer  impedance 
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For  Ztr  (Fig.  18),  the  T-parameters  of  the  two-port  network  of  tissue  properties  and 
gas  compression  can  be  derived  from  linear  circuit  analysis  (see  details  in  Appendix  A): 


P  cw 

>' 

1  +  —  z, 
Zg 

_Qcw_ 

£ 

II 

Q'. 

JL  i 

_Q'_ 

L  z* 

(31) 


It  may  be  noticed  that  Ts  (Fig.  17)  is  the  transmission  matrix  from  trachea  to  alveoli. 
The  transmission  matrix  from  the  alveoli  towards  the  trachea,  Tcw_o,  can  be  expressed  in 
terms  of  the  parameters  in  the  matrix  Ts  as 


~Po 

ts22 

ts\2 

■ Pc' 

—  Tew  _  o  • 
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ts\\ 

_Q°. 

(32) 


Therefore,  the  complete  transmission  matrix  from  alveoli  to  trachea  is 


P  cw 

'Po' 

ii 

Tew  _  0  • 

_Qcw_ 
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Combining  Eqns  3 1  and  33,  we  have 


(33) 
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(34) 


The  transfer  impedance  Ztr  is  thus  given  by 


P.= o  =(1 + — ) ' Ln  +  • /jl  1 

Zg 


(35) 


3.2.5  Pressure  ratios 

The  ratios  of  the  pressure  at  the  ends  of  airway  branches  of  a  given  generation, 
divided  by  the  applied  pressure  can  be  computed  using  the  same  scheme.  In  general,  the 
airway  opening  is  terminated  by  a  load  impedance,  Zl  (Fig  .19). 

Zi„  at  the  end  of  each  stage  is  the  load  impedance  of  the  next  stage.  For  example, 
Zin(l)  in  Fig  .19  is 

(36) 

f/>2i(l)  •  Zl  +  tpn(l) 
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where  tpll(l),  tpl2(l),  tp2,(l),  and  tp22(l),  are  parameters  in  the  transmission  matrix  Tpl 
obtained  in  section  3.2.2. 

Zin  of  the  next  stage,  Zin(2)  is 

.  .  _  tpu{2)  •  Zm(l)  +  tpn(2)  (37) 

H  '  ~  /pji(2)  •  Zm(l)  +  tpn(2) 

where  tpI1(2),  tplI(2),  tp21(2),  and  tp22(2)  are  parameters  in  the  transmission  matrix  Tp2. 

This  is  also  a  recursive  process.  Therefore,  the  input  impedances,  Zjn(i),  i=l,  2,  ...,  k 
can  be  computed  using  a  recursive  algorithm: 

2in(\\  =  +  (3  8) 

(p2i(/)  *  Zin(j  —  1)  *+*  tp2i(l) 

where  tpn(i),  tpi2(i),  tp2i(i),  and  tp22(i)  are  parameters  in  ith  transmission  matrix  Tp,  for  ith 
generation : 


cosh(y(i)  *  l(i)) 
m(i)  •  sinh(y(i) » l(i)) 


Zo(i)  •  sinh(y(i)  •  l(i)) 

m(0 

cosh(y(i)  •  l(i)) 


(39) 


where  m(i)  is  the  total  number  of  airway  segments  of  the  7th  generation,  l(i)  is  the  length, 
y(i)  is  the  propagation  coefficient,  and  Zo  is  the  characteristic  impedance. 
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Letting  Pcw=l,  the  pressure  ratio  at  each  generation  is  equal  to  the  pressure  at  the 

pfi) 

corresponding  generation,  namely,  R(i)  =  — —  =  P(i).  Multiplying  both  sides  of  Eq.  (31) 

Pew 

with  the  inverse  of  Tti,  we  have: 


m. 


1  -Zu 


.  i Zti  n 

1+ —  Q* 


where  P(k)  and  Q(k)  correspond  to  P’  and  Q’  in  Eq.  (31),  respectively.  The  pressure  P(k) 


P(k)  =  Pc»-Zti-Qc»  =  \-Z«- Qc 


where  Qcw  is  computed  as  1  divided  by  Zi„_cw  since  Pcw=l;  and  Zjn_cw  is 


(1  +  —)-Zin(k)  +  Zti 

MV ~ 


In  general,  the  pressure  at  generation  /  can  be  calculated  using  the  recursive 
algorithm: 


P(i)  =  tP  1 1  (/  + 1)  •  P(i  + 1)  -  TPn(i  + 1) 


PQ+ 1) 

Z„(i+ 1) 


Load  impedance  Zl  is  zero  when  the  mouth  is  open  and  approaches  infinity  when  the 


mouth  is  closed. 
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3.3  Input  and  transfer  impedance  for  asymmetrically  branching  tree 


structure 

3.3.1  Self-similarity  algorithm  for  calculating  input  impedance  Zm 

Zinofthe  asymmetrical  model  is  computed  using  self-similarity  algorithms  [1,  15].  It 
starts  with  the  most  peripheral  airway  and  marches  up  the  longest  pathway  to  the  airway 
opening  using, 


Zt  +  Zo  tanh(yl(n)) 

1  +  (Zt  /  Zo)  tanh(yl(n)) 

where  ZT  is  the  parallel  combination  of  the  two  daughter  branches  of  order  n. 


(44) 


Zt  =  Zin(n  -  1)//Zin(n  -  1  -  A(n)) 


(45) 


When  n=l  (terminal  airway),  Zt  is  the  load  impedance  of  the  series  tissue  impedance 
in  parallel  with  alveolar  gas  compression  capacitance  as  described  in  Eq.  (28). 


3.3.2  Principle  of  superposition 

Fluid  mechanics  in  the  airways  due  to  small  amplitude  oscillations  are  considered  to 
behave  linearly  and  thus  obey  the  superposition  theorem.  The  superposition  theorem  says 
that  [2]: 

In  any  linear  network  containing  several  sources,  the  voltage 
(analogous  to  pressure)  across  or  the  current  (analogous  to  flow)  through 
any  branch  can  be  calculated  by  adding  algebraically  all  the  individual 
voltages  or  currents  caused  by  each  independent  source  acting  alone,  with 
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all  other  independent  voltage  sources  replaced  by  short  circuits  and  all  the 
other  independent  current  sources  replaced  by  an  open  circuit. 

Let  us  first  consider  applying  the  pressure  sources  to  the  terminal  airways  of  the 
simple  tree  network  in  Fig.  20.  The  flow  due  to  a  pressure  applied  to  the  1st  terminal 
airway  alone  (with  pressure  sources  at  2nd  and  3rd  terminals  grounded)  is  Flowj.  The  flow 
due  to  a  pressure  applied  to  the  2nd  terminal  alone  is  Flow2,  and  the  flow  due  to  a  pressure 
applied  to  3rd  terminal  alone  is  Flow3.  According  to  the  superposition  theorem,  the  total 
flow  in  each  of  the  airway  segments  due  to  pressure  sources  simultaneously  applied  to  all 
of  the  terminal  airways  is  simply  the  sum  of  Flow  ^  Flow2,  and  Flows,  Ztr  of  the  whole 
network,  where  all  pressure  sources,  P,  are  identical,  is: 

^'r  =  ZFlow*  (46) 


43 


The  superposition  theorem  enables  us  to  compute  Ztr  by  considering  each  single 
pressure  source  one  by  one.  That  is,  it  is  not  necessary  to  consider  the  network  connection 
globally  as  we  did  to  analyze  Ztr  in  the  symmetrical  tree  structure. 

If  the  input  impedances  looking  from  bottom  up  along  an  airway  path  are  available,  it 
is  quite  straightforward  to  compute  flows  entering  all  the  airway  segments  on  this  path. 
The  output  pressure  of  this  segment  can  be  calculated  from  the  two-port  network 
parameters,  which  are  the  input  pressure  of  its  parent.  Therefore,  marching  up  the  path 
with  this  recursive  process,  we  can  obtain  the  flow  due  to  the  pressure  applied  at  the 
bottom  of  this  path.  The  remaining  question  is,  how  do  we  calculate  input  impedances 
looking  up  from  the  end  of  all  the  airway  segments? 

3.3.3  Analysis  of  the  structure  of  Horsefield’s  model 

Generally,  from  each  different  terminal  airway,  a  different  tree  structure  is  seen.  In 
order  to  calculate  the  input  impedance  looking  up  from  that  terminal  airway,  the 
corresponding  tree  structure  should  be  used.  Let  us  consider  the  Horsefield  model  for  the 
human  lung.  The  tree  is  constructed  based  on  the  following  rule:  the  n,h  generation,  the 
“parent”,  bifurcates  into  two  daughters;  the  left  daughter  is  always  n-1,  and  the  right 
one  is  n-l-A(n)  where  A(n)  is  the  recursion  index. 

For  the  sake  of  simplicity,  let  us  consider  the  shortest  pathway  in  the  human  lung 
model  (Fig.  21).  The  impedance  looking  up  the  tree  towards  the  airway  opening  of  the  right¬ 
most  path  is  equivalent  to  the  input  impedance  of  the  network  shown  in  Fig.  22. 
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Note  that  Zink  is  the  Zin  looking  into  the  k,h  generation  towards  the  alveoli,  all  of 
which  have  been  calculated  and  stored  when  computing  Zin. 


Z3 6 


Figure  22  Conceptural  network  model  for  shortest  pathway  in  Horsefield’s 
model  of  the  human  lung 

The  input  impedance  of  its  adjacent  path  is  equavalent  to  the  the  input  impedance  of 
the  circuit  network  shown  in  Fig.  23.  Note  that  the  shaded  parts  in  Fig.  22  and  Fig.  23  are 
the  same,  which  indicates  that  we  do  not  have  to  re-calculate  the  impedance  along  the  entire 
path.  In  fact,  if  the  impedance  looking  into  A  is  stored,  only  the  combination  of  five- 
elements  {Zn,  Zinn,  Z13,  Zina,  ZinA}  instead  of  seventeen-element  path  needs  to  be 
computed  in  this  case. 
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Figure  23.  Conceptual  network  model  for  the  adjacent  pathway  of  the  shortest  in  the 
Horsefield  model 

3.3.4.  Horizontal-decomposition  of  the  tree  structure 

—  symmetrical  zone  and  asymmetrical  zone 
The  Horsefield  model  describes  the  branching  pattern  and  dimensions  of  all  the 
airways  in  the  lung  down  to  and  including  the  respiratory  bronchioles.  For  the  sake  of 
computational  convenience,  it  can  be  decomposed  into  two  structural  zones  (Fig.  24).  In  the 
asymmetrical  zone  from  order  36  to  order  12  the  branching  is  asymmetrical,  all  of  which 
have  two  different  daughter  branches  (except  order  12).  The  reason  to  include  order  12  in 
the  asymmetrical  zone  is  that  order  14  has  a  recursion  index  of  2  and  has  two  daughter 
branches  of  order  13  and  order  1 1,  respectively,  such  that  order  12  is  “skipped”  in  this  case. 
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In  addition,  order  13  has  a  recursion  index  of  1,  which  results  in  daughter  branches  of  order 
12  and  11.  The  tree  structure  becomes  completely  symmetrical  from  order  11  to  order  1.  In 
fact,  order  14  to  order  12  can  be  considered  as  the  transition  zone.  Since  there  is  no 
computational  benefit  to  having  this  extra  level  of  decomposition,  it  is  combined  with  the 
asymmetrical  zone.  As  one  can  see,  in  the  symmetrical  zone  (order  1 1  to  order  1),  the 
structure  is  symmetrical,  which  has  the  computational  simplicity  and  beauty  as  discussed  in 
section  3.2. 

3.3.5.  Vertical-Decomposition  of  the  asymmetrical  part 

In  order  to  solve  the  problem  of  limited  memory  versus  large  storage  needs,  further 
decomposition  is  required,  namely,  vertical  decomposition  (Fig.  25). 


Figure  25.  Vertical  decomposition  of  the  tree  structure 


50 


When  we  compute  the  input  impedances  looking  from  bottom  up  at  most  nodes  in 
subtree  1,  the  detailed  structural  information  of  subtree  2  and  subtree  3  are  not  relevant. 
Therefore,  we  would  rather  not  use  memory  to  store  any  parameters  related  to  them.  This  is 
more  critical  when  the  computational  resources  are  limited  while  the  tree  contains  a  larger 
number  of  segments.  The  three  subtrees  are  almost  disjoint  in  that  subtree  1  and  subtree  2 
only  have  order  35  and  order  36  in  common,  and  subtree  1  and  subtree  3  only  have  order  36 
in  common.  Hence  there  is  little  repetition  of  computation  and  storage. 

3.3.6.  Unique  indexing  binary  tree  representation  of  Horsefield  model 

In  the  unique-index  binary  tree  (Fig.26),  the  indices  start  from  1,  which  is  the  index  for 
the  root  of  the  tree.  If  the  index  for  parent  branch  is  n,  the  indices  of  the  two  daughter 
branches  (if  they  exist)  will  be: 

(47) 

(48) 


Index  of  left  daughter  =  2n 
Index_of_right_daughter  -  2n+I 


51 


One  of  the  advantages  of  unique-index  tree  representation  is  that  all  airway  branches 
have  their  own  indices  such  that  all  of  them  can  be  uniquely  identified  and  computed 
individually.  Another  more  important  advantage  is  that  the  numbering  itself  “indicates”  the 
position  of  the  branch  in  the  entire  network  in  that  its  parent,  grandparent,  up  to  the  root  of 
the  tree  can  be  derived  from  its  index.  The  parity  of  the  index  indicates  whether  it  is  on  the 
left  or  on  the  right,  and  therefore  its  neighbor  can  be  identified  as  well.  This  is  illustrated  in 
Fig.  27.  For  a  given  branch  whose  index  is  13,  let  us  find  the  path  to  get  here: 

1)  Since  13  is  odd,  this  branch  must  be  a  right  daughter  branch,  and  its  neighbor 

13 

(on  the  left)  is  12.  The  index  of  its  parent  is  —  ("[xj:  “floor  of  x”,  the  greatest 
integer  <  x),  which  is  6. 

6 

2)  6  is  even,  so  it  is  on  the  left  and  its  right  neighbor  is  7.  Its  parent  is  —  =3. 

3)  3  is  odd,  so  it  is  on  the  right  and  its  left  neighbor  is  2.  Its  parent  is  1,  which  is  the 
root  of  the  tree. 

Note  that  if  the  index  /  of  the  branch  is  odd,  its  neighbor’s  index  is  i  -1.  Otherwise  its 
neighbor’s  index  is  i  +  1; 
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This  turns  out  to  be  a  powerful  feature  when  implementing  the  algorithm  for 
computing  the  transfer  impedance  Ztr  since  searching  is  straightforward. 

However,  the  actual  morphometric  data  are  associated  with  Horsefield  order  number. 
Therefore,  it  is  necessary  to  map  the  Horsefield’s  model  onto  unique-indexing  binary  tree 
representation.  For  a  given  unique  index,  the  corresponding  Horsefield  order  number  will 
indicate  where  to  find  the  morphometric  data.  The  two  indices  are  mapped  in  Table  9: 


Table  9.  Mapping  Horsefield  model  to  Unique  index  representation 


Unique  index  representation 

Horsefield  model 

Root  of  the  tree 

1 

N  (total  number  of  orders) 

For  a  given  branch 

n 

k 

Left  daughter 

2n 

k-1 

Right  daughter 

2n+l 

k-l-A(k) 

3.3.7.  Recursive  algorithm  for  searching  all  the  airway  paths 

For  a  given  destination,  i.e.  Kth  order,  the  following  algorithm  is  used  to  obtain  all  the 
possible  pathways  when  starting  from  the  nth  order. 

For  left  branch  (node): 

Function  lnod((n,  k,  path ,  flag) 

if(n>K)  %  if  haven’t  reached  the  destination  yet 

n=n-l/  %  keep  going  down  left  further 

if  (n==K)  %  if  reached  the  destination ,  record 
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print  [path,  n] 
print  [flag,  0] 
elseif  (n>K)  %  keep  going 

lnod  (n,K,  [path,n]  ,  [flag ,0] )  ; 
mod  (n,K,  [path,n]  ,  [flag,0])  ; 
end; 
end; 


For  right  branch  (node): 

Function  mod  (n ,  k,  path,  flag) 

if(n>k)  %  if  haven't  reached  the  destination  yet 

n=n-l-delta  (n)  ;  %  keep  going  down  right  further 
if(n~k) 

print  [path,  n] 
print  [ flag,  1] 
elseif  (n>k) 

lnod  (n,k,  [path,n] , [flag,l] )  ; 
mod  (n,k,  [path,n]  ,  [flag,l] )  ; 
end; 

end; 

The  algorithm  is  fairly  simple.  The  speed  of  the  algorithm  depends  on  the  maximum 
pathway  length  (how  many  to  go  through)  and  the  recursion  indices. 
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3.3.8.  Transfer  impedance 

In  the  asymmetrical  zone  of  the  Horsefield’s  model,  all  of  the  paths  end  with  order 
12, 13  or  14  (Fig.  28).  The  number  of  paths  and  the  number  of  nodes,  which  is  equal  to  the 
number  of  different  input  impedances,  of  sub-tree  1  are  listed  in  Table  10.  Starting  from  the 
top  of  each  path  we  compute  the  input  impedances  looking  from  bottom  up.  If  the  input 
impedance  of  the  node  (the  end  of  the  corresponding  airway  branch)  has  already  been 
calculated,  we  skip  it  and  look  at  the  next  one  on  the  path  until  the  end  is  reached. 


Figure  28.  Illustration  of  the  three  possible  ways  which  all  the  paths  in  three  subtrees 
can  end  with. 
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Table  10.  The  number  of  paths  and  nodes  in  Sub-tree  1 


Path  ends  at: 

Order  12 

Order  13 

Order  14 

Number  of  paths 

907 

476 

345 

Total  number  of  paths 

1728 

Total  number  of  nodes 

2633 

As  one  can  see  that  all  these  paths  have  many  nodes  in  common.  We  avoid  repeating 
the  computation  of  the  input  impedance  of  the  same  node  by  determining  whether  the 

unique  index  of  this  node  has  been  marked. 

For  the  symmetrical  zone  (from  order  1 1  to  the  most  peripheral  terminal  airways),  the 
methods  discussed  in  section  3.2  are  applied  to  obtain  the  equivalent  two-port  network 
block  Tsym  combining  airways  from  order  11  to  1  and  the  tissue  properties  (Fig.  29).  The 
impedance  looking  up  from  the  last  node  of  a  given  path,  in  parallel  with  the  input 
impedance  of  the  neighbor  of  order  1 1th,  serves  as  the  load  impedance  of  this  two-port 
network. 
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Figure  29.  Two-port  network  representation  of  symmetrical  zone 

The  output  pressure  of  the  two-port  network,  TSym,  is  the  equivalent  pressure  source 
to  the  corresponding  path  in  asymmetrical  zone.  Here,  m  in  Eq.  (22)  is  2n  n,  where  n  is  the 
order  number  in  Horsefield’s  model.  Nends  in  Eq.  (30)  is  210,  instead  of  total  number  of 
terminal  airways. 

Therefore,  proceeding  up  the  path  in  the  asymmetrical  zone,  the  pressure  at  the  end 
of  trachea  and  the  flow  at  the  airway  opening  can  be  computed.  The  sum  of  all  flows  due  to 
all  of  the  identical  applied  pressure  sources  is  the  total  flow  in  the  trachea,  EFlowk. 
Dividing  the  pressure  by  this  flow  sum  gives  the  transfer  impedance, 


7  =  — - — 

tr  XFlowfc 


(49) 
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3.3.9.  Pressure  ratio  distribution  along  the  shortest  path 

To  compute  the  distribution  of  pressure  ratios  along  the  shortest  path  in  the 
asymmetrical  zone,  a  tractable  organization  of  the  entire  tree  network  is  needed.  The  entire 
tree  is  re-organized  (pruned)  with  respect  to  the  shortest  path  (Fig  30)  into  seven  subtrees. 
The  equivalent  flow  sources  from  the  corresponding  sub-trees  are  identified  as  flow_l 
through  flow_7.  The  nodes  where  the  pressure  ratios  are  computed  are  labeled  accordingly 
(1  through  7).  The  flow  from  a  given  source  (i.e.,  flow_n)  at  each  (kth)  of  seven  nodes  will 
be  downward  if  n  <  k  and  upward  if  n  >  k.  The  flow  due  to  all  the  flow  sources  is  then 
computed  using  the  superposition  theorem.  Since  the  impedance  at  each  node  is  known,  the 
pressure  at  that  node  can  be  computed  in  a  straightforward  way  and  normalized  to  the 
applied  pressure. 
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4.  Results  and  Discussion 


4.1  Symmetrical  model  versus  asymmetrical  model 

4.1.1  Input  Impedance  Zi„  of  the  human  lung  models 

Weibel’s  symmetrical  model  (Table  1)  and  Horsefield’s  asymmetrical  model  (Table3) 
both  combined  with  non-rigid  airway  wall  parameters  (Table  5  and  Table  4,  respectively) 
were  used  to  simulate  Zin  (Fig.  31).  Note  that  neither  model  included  the  upper  airways. 
The  first  anti-resonance  in  Weibel’s  symmetrical  airway  network  occurs  at  a  lower 
frequency  than  it  does  in  Horsefield’s  asymmetrical  airway  network.  The  dominant  factor 
that  accounts  for  the  differences  in  the  frequency  of  this  anti-resonance  is  the  length 
difference  in  the  trachea  rather  than  the  symmetry  of  the  airway  network.  In  Weibel’s 
model,  the  trachea  is  12  cm  in  length  while  it  is  10  cm  long  in  Horsefield’s  model.  This 
complies  with  the  fact  that  the  frequency  of  the  acoustic  resonances  and  anti-resonances  in 
a  tube  is  inversely  related  to  its  length.  The  experimental  Zin  data  published  by  Habib  et  al. 
[1]  are  shown  in  Fig.  31  (b)  for  comparison  to  the  model  predicted  Zm.  The  Z-m  predictions 
from  both  symmetrical  and  asymmetrical  airway  network  compare  closely  to  the  measured 
data. 
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Figure  31.  (a)  Z-m  based  on  Weibel’s  symmetrical  airway  network  (*)  and  Z-m  based  on 
Horsefield’s  asymmetrical  airway  network  (o) 

(b)  Subglottal  (i.e.,  without  upper  airways)  Zm  data  (•)  measured  in  adult 
human  (Habib  et  al.  [1]) 


4.1.2  Total  respiratory  input  impedance  including  upper  airways 

The  total  respiratory  Zin  including  upper  airways  was  simulated  and  compared  to 
measured  data  [22]  (Fig.  32).  The  resonant  peak  is  shifted  to  lower  frequencies  compared  to 
subglottal  Zin  due  to  the  added  airway  length.  Summing  the  lengths  of  the  upper  airways 
from  Table  7,  we  have  a  total  length  of  23.2  cm  for  the  upper  airways. 


(a)  (b) 


FREQUENCY 

Real  (top)  and  imaginary  (bottom)  parts  of  human  input 
impedance  averaged  from  all  8  subjects  (Table  1).  Mean  (*)  ±  SD 
( - ),  Units  of  impedance  are  cmHjO*!-1^**1. 


Figure  32.  (a)  Total  respiratory  Zm.  (b)  Zm  data  measured  in  adult  human  (Dorkin,  H.  L.  et 


al.  [22]) 
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4.1.3  Ztr  from  WeibePs  model  versus  that  from  Horsefield’s  model. 

Ztr  is  usually  measured  with  the  setup  shown  in  Fig.  2  (a).  The  forced  oscillation  is 
placed  around  the  chest  wall,  and  flow  is  measured  at  the  airway  opening.  Since  the  mouth 
of  a  subject  is  open,  the  impedance  to  flow  at  the  airway  opening  is  extremely  small  and 
modeled  as  a  short  circuit  (Fig.  18).  The  upper  airways  influence  the  measurement,  and 
therefore,  they  were  included  in  the  simulations.  Ztr  was  computed  from  1  up  to  72  Hz  in 
both  the  symmetrical  and  asymmetrical  models  and  compared  with  measured  data 
published  by  Lutchen  et  al  [7]  (Fig.  33). 

To  explore  Ztr  at  higher  frequencies  in  more  detail,  we  simulated  it  up  to  1000  Hz 
using  the  symmetrical  and  asymmetrical  models  (Fig.  34).  They  are  very  similar  except  in 
the  symmetrical  airway  model  the  anti-resonance  occurs  at  a  lower  frequency  than  in  the 
asymmetrical  airway  model.  Again,  this  is  most  likely  because  of  the  longer  trachea  in  the 
symmetrical  model. 
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Figure  33.  Ztr  of  the  symmetrical  human  lung  model  (•),  Ztr  of  the  asymmetrical 
human  lung  model  (o),  and  measured  data  (average  ‘  ’  ±  SD  at  low  frequencies.. 
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Figure  34  Ztr  of  the  Weibel’s  symmetrical  human  lung  model  (•)  and  Ztr  of  the 
Horsefield’s  asymmetrical  human  lung  model  (-)  at  high  frequencies 
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4.1.4  Pressure  ratio  distribution  with  forced  oscillation  imposed  at  the  airway  opening 
It  has  been  reported  that  in  Horsefield’s  dog  lung  model  the  serial  pressure 
distributions  are  frequency-dependent  and  pressures  in  the  peripheral  airways  and  alveoli  at 
some  frequencies  can  be  larger  than  at  the  airway  opening  [20],  We  were  interested  to 
know  whether  this  was  also  the  case  in  the  asymmetrical  human  lung  model.  First,  we 
simulated  the  pressure  ratio  at  low  frequencies  from  1  to  100  Hz  along  the  longest  pathway 
in  the  human  lung  model  (Fig.  35  (a)).  We  found  that  the  ratio  of  alveolar  pressure,  Paiv,  to 
airway  opening  pressure,  Pao,  is  greater  than  unity  below  5  Hz,  which  agrees  with  the 
findings  of  Jackson  et  al  in  the  dog  lung  [20],  The  maximum  of  the  pressure  ratio  is  about 
1.27  and  occurs  at  4  Hz.  Pressures  were  amplified  below  5  Hz  in  all  the  airway  orders  from 
the  most  peripheral  airway  to  the  upper  airways.  The  maxima  of  the  pressure  ratios  in  all 
airway  orders  occur  at  4  Hz  (Fig  35  a),  and  at  this  frequency  the  pressure  ratio  reaches  its 
maximum  in  airway  order  15  (Fig  35  b).  To  compare  our  predictions  in  the  human  model  to 
those  of  Jackson  et  al  (20),  we  predicted  the  pressure  ratios  in  the  trachea  and  alveoli  (Fig 
35  c,  d). 
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Figure  35.  Pressure  ratio  distribution  at  low  frequencies,  (a)  Pressure  ratio  distribution 
versus  frequency  and  airway  order,  (b)  Pressure  ratio  versus  airway  order  at  4  Hz,  where  the 
maximum  occurs,  (c)  Pressure  ratio  in  the  trachea  (36th  order),  (d)  Pressure  ratio  in  the 


alveoli  (1st  order) 


We  further  explored  the  distribution  of  the  pressure  ratios  at  frequencies  up  to  2200  Hz 
(Fig.  36  (a))  and  found  that  pressures  in  the  central  airways  can  be  6  to  10  times  higher  than 
Pao  at  certain  frequencies,  i.e.,  approximately  400  Hz,  800  Hz,  and  1600Hz.  Fig.  36  (b),  (c) 
and  (d)  illustrate  in  detail  the  pressure  ratios  in  the  upper  airway  (41st  order),  the  trachea 
(36th  order)  and  the  peripheral  airway  (1st  order),  correspondingly.  Quite  different  from 
what  was  found  at  low  frequencies,  the  pressure  ratios  are  much  higher  in  the  central 
airways  than  in  the  peripheral  airways.  There  is  no  pressure  amplification  in  the  peripheral 
airways  for  f  >  5  Hz. 
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Figure  36.  Pressure  ratio  distribution  at  high  frequencies,  (a)  Pressure  ratio  distribution 
versus  frequency  and  airway  order,  (b)  Pressure  ratio  in  the  upper  airway  (41st  order),  (c)  in 
the  trachea  (36th  order),  (d)  in  alveoli  (1st  order) 


4.1.5  Pressure  ratio  distribution  with  forced  oscillation  imposed  at  the  chest  wall 

To  investigate  the  potential  influence  of  low  frequency  sound  (LFS)  on  divers,  we 
simulated  the  pressure  ratios  versus  frequency  and  airway  order.  The  pressure  distribution 
was  simulated  for  a  boundary  condition  where  Z=oo  at  the  airway  opening,  i.e.,  closed 
mouth  and  closed  nasal  passage.  This  is  analogous  to  having  an  open  circuit  at  the  airway 
opening.  The  pressure  ratios  as  a  function  of  frequency  (up  to  250  Hz)  and  airway  order 
were  first  simulated  using  the  symmetrical  model  and  no  pressure  amplification  was  found 
(Fig.  37  (d)).  In  fact,  in  all  airways  the  pressure  decreases  significantly  as  frequency 
increases.  The  maximum  pressure  ratio  occurs  at  the  lowest  frequency  (used  1  Hz)  and  the 
pressure  ratio  is  relatively  independent  of  airway  generation.  Intuitively,  it  is  difficult  to 
understand  why  there  would  be  such  a  small  pressure  change  along  the  airways.  Since 
pressure  is  related  to  impedance  and  flow,  it  is  natural  to  study  impedance  and  flow 
distribution  versus  airway  generation.  In  the  equivalent  series  two-port  network  (Fig.  19), 
the  impedance  (Zcw.ao)  looking  from  the  chest  wall  up  the  airway  tree  as  a  function  of 
airway  order  at  1  Hz  is  shown  in  Fig.  37  (a).  The  flow  (Fig.  37  b)  as  a  function  of  airway 
order  decreases  from  the  alveoli  to  the  airway  opening  while  Zcw-ao  increases.  This  is  a 
reciprocal  process  resulting  in  little  change  in  the  pressure  ratio  (Tig.  37  c). 
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The  pressure  distribution  in  the  asymmetrical  airway  model  depends  not  only  on  the 
airway  order,  but  also  on  the  path  taken  to  a  particular  airway  segment.  Simulating  the 
pressure  ratio  distribution  over  the  entire  airway  network  using  the  superposition  model  is 
not  a  trivial  task.  Thus  for  simplicity,  we  took  the  shortest  pathway  as  the  representative 
pathway  and  simulated  the  pressure  ratios  along  this  path.  To  further  simplify  the 
simulation,  the  airways  below  the  14th  order  were  not  included  since  these  airways  branch 
symmetrically  and  nothing  interesting  was  found  in  peripheral  airways  in  the  symmetrical 
airway  network  (Fig.  38  (d».  The  pressure  ratio  is  greater  than  unity  (i.e.,  pressure 
amplification)  at  ~40Hz  in  the  trachea  and  upper  airways  (Fig.  38  (d)).  We  found  a 
significant  drop  in  the  pressure  ratio  going  from  the  trachea  into  the  next  airway  order.  To 
study  this  further,  we  plotted  the  input  impedances  looking  from  the  chest  wall  towards  the 
airway  opening  (Zcw_ao)  and  flow  as  a  function  of  airway  order.  A  large  impedance  drop 
was  found  between  the  trachea  (36th  order)  and  the  34th  order  airway.  The  difference 
between  the  symmetrical  and  asymmetrical  models  is  that  the  magnitudes  of  both  Zcw-ao  and 
flow  are  increasing  when  proceeding  from  the  peripheral  airways  to  the  upper  airways.  This 
may  explain  why  the  pressure  ratio  distribution  in  the  asymmetrical  model  looks  very 
different  from  that  of  the  symmetrical  model. 
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Figure  38.  Pressure  ratio  distribution  in  the  Horsefield’s  asymmetrical  airway  network, 
(a),  (b)  and  (c)  Zcw-ao ,  flow,  and  pressure  versus  airway  order  at  1  Hz. 

(d)  pressure  ratio  distribution  as  a  function  of  frequency  and  airway  order. 
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4.2  Strain,  wave  propagation  velocity  and  wall  admittance  distributions  in 
asymmetrical  airway  network 

4.2.1  Strain  as  a  function  of  frequency  and  airway  order 

It  is  not  clear  whether  barotrama  in  the  airways  is  due  to  excessive  pressure  (stress) 
or  deformation  (strain)  of  the  airway  wall  (i.e.,  the  motion  or  displacement  of  the  airway 
wall).  For  example  consider  two  different  airways  that  are  exposed  to  a  transmural  pressure 
oscillation  of  the  same  magnitude;  if  one  has  a  smaller  magnitude  of  impedance  it  would 
experience  a  larger  strain.  Strain  in  a  tissue  is  usually  defined  as 


e=(L  -  Lo)/L0 


(50) 


where  Lo  is  the  tissue’s  resting  length  and  L  is  its  deformed  length.  However,  here  we  will 
consider  the  volume  strain  of  an  airway  segment  where  volume  strain,  €v  is  defined  as 


e  v  =  (V-Vo)/Vo 


(51) 


where  Vo  is  the  resting  volume  of  the  airway  segment  and  V  is  its  maximally  deformed 
volume.  The  magnitude  of  the  time  derivative  of  an  airway  segments  ev  is  the  magnitude  of 
its  volume  flow  and  is  related  to  its  wall  impedance  and  the  transmural  pressure  by, 


dev/dt  =  Ptm/Zaw. 


(52) 
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The  pressure  and  thus  the  flow  will  be  a  sine  function,  or 


dev/dt  =  Asin(©  t).  (53) 

The  magnitude  of  the  segment’s  volume  change  AV,  is 

AV  =  J  flow  dt  =  J  A-sin(©t)dt  =  A-cos(©t)  /  ©  (54) 

Thus  the  amplitude  of  ev  is  given  by 

€v=  A/©  (55) 

The  values  of  ev  within  the  airway  tree  depend  on  the  amplitude  of  the  applied  pressure  (i.e., 
Pao  when  applied  at  the  airway  opening,  or  PbS  when  applied  at  the  body  surface).  Thus,  we 
predicted  €v  for  a  unity  Pao  and  Pbs  (Fig.  40).  The  strain  is  normalized  to  its  maximum  value 
and  therefore  the  maximum  strain  is  unity.  The  volume  strain  ev  is  a  maximum  at 
approximately  30  Hz  and  occurs  in  the  central  airways. 
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4.2.2  Wave  propagation  velocity  as  a  function  of  frequency  and  airway  order 


In  rigid-walled  tubes,  the  velocity  of  wave  propagation  is  equal  to  the  free-field  speed 
of  sound  (34,318  cm/s  in  room  air).  However,  in  nonrigid-walled  tubes,  wave  propagation 
velocity  becomes  frequency  dependent  [21].  The  wave  propagation  velocity  distribution  as 
a  function  of  frequency  and  airway  order  are  shown  in  Fig.  41. 


order  Frequency  (Hz) 


Figure  41.  Wave  propagation  velocity  distribution  as  a  function  of  frequency  and  airway 
order. 
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4.2.3  Relationship  between  strain  and  wave  propagation  velocity 

If  the  shunt  admittances  of  a  nonrigid  walled  tube  are  significantly  under-damped 
(i.e.,  their  resistance  is  small),  sound  may  not  be  propagated  at  or  near  its  resonant 
frequency  because  wave  propagation  velocity  could  become  less  than  zero  (Fig.  42).  At 
frequencies  where  c<0,  there  would  be  extremely  large  volume  strains  and  no  energy  would 
be  propagated  since  it  is  all  dissipated  in  wall  motion. 


Figure  42.  Wave  propagation  velocity  with  no  propagation  at  certain  frequencies 

To  explore  the  relationship  between  strain  and  wave  propagation  velocity,  we  plotted 
them  together  in  Fig.  43.  These  are  the  projections  on  the  frequency  axis  of  their  3-D  plots, 
shown  in  Fig.  40  and  Fig.  41,  respectively.  Different  curves  represent  different  airway 
orders.  As  we  expected,  the  regions  of  maximum  stain  correspond  to  the  regions  of 
minimum  wave  propagation  velocity,  as  indicated  by  dashed  line  boundaries. 
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Figure  43.  The  relationship  between  the  strain  and  wave  propagation  velocity. 
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4.2.4  Relationship  between  strain  and  airway  wall  admittance 

We  took  the  airway  segment  (order  41st)  in  which  strain  achieves  its  maximum,  as  the 
representative  to  study  the  relationship  between  strain  and  airway  admittance.  In  Fig.  44  are 
plotted  the  strain  as  a  function  of  frequency  (a),  the  magnitude  of  wall  admittance  per  unit 
length  Yw  as  a  function  of  frequency  (b),  the  real  part  of  Yw  (c),  and  the  imaginary  part  of 
Yw  (d).  The  maximum  in  strain  and  magnitude  of  Yw  occur  at  approximately  the  same 
frequencies,  and  the  wall  admittance  shows  an  anti-resonance. 
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Figure  44.  Relationship  between  strain  and  airway  wall  admittance  (order  41  is  taken  as  the 
representative.) 


4.3  Simulations  in  asthmatic  lung  with  comparison  to  those  in  normal  lung 

The  Horsefield’s  asymmetrical  model  with  upper  airways  was  used  to  simulate  Zin 
and  Ztr  in  lungs  with  homogeneous  peripheral  airway  obstruction  and/or  increase  peripheral 
airway  wall  thickness  (a  model  of  asthma). 

4.3.1  Comparison  of  Zin  and  Ztr  of  normal  lung  and  with  asthma 

In  the  peripheral  airways  (i.e.,  d  <  2  mm),  we  decreased  diameters  and/or  increased 
wall  thickness.  Two  different  degrees  of  constriction  were  simulated,  40%  and  60% 
decrease  in  airway  diameters  and  two  different  increases  in  wall  thickness  were  simulated, 
40%  and  60%.  When  the  airways  become  constricted  the  frequency  at  which  the  acoustic 
anti-resonance  occurred  (at  -180  Hz)  did  not  change,  but  the  peak  in  the  real  and  imaginary 
parts  decreased  as  the  level  of  bronchoconstriction  increased  (Fig.  45).  The  real  part  also 
became  negatively  frequency  dependent  at  low  frequencies  with  the  highest  level  of 


bronchoconstriction. 


Figure  45.  Comparison  of  Zin  of  normal  lung  (-),  asthmatic  lung  with  40%  constriction 
and  asthmatic  lung  with  60%  constriction  (•). 
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We  also  simulated  Ztr  in  healthy  and  asthmatic  lung  (Fig.  46)  with  40%  or  60% 
constriction,  and/or  40%  or  60%  increase  in  airway  wall  thickness.  The  results  of  normal 
lung  are  compared  in  Fig.  46.  From  the  magnitude  of  Ztr,  airway  constriction,  compared  to 
normal  lung,  increases  the  transfer  impedance.  When  the  airway  wall  thickness  was 
increased  by  40%  or  60%  in  conjunction  with  airway  constriction,  there  was  little 
difference  in  Zin  with  constriction  alone  (Fig.  46). 
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Figure  46.  Comparison  of  Ztr  of  normal  lung  (-),  asthmatic  lung  with  40%  constriction 
(— ),  and  with  40%  decrease  in  diameter  alone  (•). 
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4.3.2  Pressure  ratio  distribution  with  forced  oscillation  imposed  at  chest  wall 

We  also  predicted  the  pressure  distribution  within  the  lungs  where  the  peripheral 
airways  were  constricted  and/or  their  walls  were  thickened.  The  pressure  ratio  as  a  function 
of  frequency  and  airway  order  was  simulated  along  the  shortest  pathway  (Fig.  47). 
Compared  to  the  pressure  ratios  of  the  normal  lung  (Fig.  38),  there  is  no  significant  change 
with  constriction  and/or  increased  wall  thickness.  The  maximum  pressure  ratio  is  about  1.3, 
which  is  slightly  smaller  than  the  maximum  pressure  ratio  in  normal  lung  (about  1.6). 


Frequency  (Hz) 


Figure  47.  Pressure  ratio  distribution  versus  frequency  and  airway  order  along  the  shortest 
path  in  asthmatic  lung 
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4.4  Simulation  for  Yeh’s  model  of  the  rat  lung 
4.4.1  Zin  and  Ztr 

We  simulated  Zin  and  Ztr  for  Yeh’s  model  of  the  rat  lung.  Since  Yeh’s  model  is 
symmetrical,  we  used  the  analysis  discussed  in  section  3.2.  The  lung  volume  was  assumed 
to  bel3.7  ml  [20],  The  density,  Young’s  modulus,  and  viscosity  of  the  cartilage  and  soft 
tissue  were  assumed  to  be  the  same  as  those  in  the  human  model.  Gas  type  was  100% 
humidified  room  air  plus  CO2. 

The  results  are  shown  in  Fig.  48.  The  first  resonant  frequency  was  found  to  be 
approximately  90  Hz,  and  there  were  well-damped  anti-resonances  at  approximately  5  kHz 
and  1 1  kHz.  Both  resonant  and  anti-resonant  frequencies  occurred  at  higher  frequencies 
than  in  the  human  lung  model.  The  reason  is  that  the  airway  length  of  the  rat  lung  is  smaller 
than  that  of  the  human  lung,  which  results  in  a  shift  of  resonant  and  anti-resonant 
frequencies  to  higher  frequencies.  Similarly,  the  resonant  and  anti-resonant  frequencies  of 
Ztr  shift  to  higher  frequencies  for  the  same  reason  (Fig.  49).  The  magnitude  of  transfer 
impedance  increases  almost  monotonically  with  frequency. 
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Figure  48.  Zin  of  the  Yeh’s  symmetrical  model  of  the  rat  lung  (-)  and  measured  data 
(•)  (Jackson  et  al.  [23]) 
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4.4.2  Pressure  ratio  with  forced  oscillation  imposed  at  airway  opening 

The  pressure  ratio  distribution  with  forced  oscillation  imposed  at  airway  opening  was 
simulated  in  a  similar  manner  as  was  done  for  the  human  lung  model.  The  result  is  shown 
in  Fig.  50.  As  we  can  see  from  Fig.  50,  the  alveolar  pressure  is  higher  than  that  at  the 
airway  opening,  which  is  similar  to  the  findings  in  the  human  lung  model.  However,  the 
frequency  range  over  which  pressure  amplification  occurs  is  larger  (1  <  f  <  100  Hz) 
compared  that  in  the  human  lung. 
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Frequency  (Hz) 


Figure  50.  Pressure  ratio  distribution  with  forced  oscillation  imposed  at  the  airway 


opening 


4.4.3  Pressure  ratio  with  forced  oscillation  imposed  at  chest  wall 


The  pressure  ratio  distribution  with  forced  oscillations  imposed  at  the  chest  wall  was 
also  investigated  in  the  rat  lung  model  using  the  methods  discussed  in  section  3.2  (Fig.  51). 
The  pressure  ratio  is  similar  to  that  in  the  symmetrical  human  lung  model.  This  might  be 
due  to  the  symmetry  property  of  both  airway  networks. 


Figure  51.  Pressure  ratio  distribution  in  the  Yeh’s  rat  lung  model  with  forced 
oscillation  imposed  at  chest  wall. 
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5.  Summary 

We  predicted  input  (Zin)  and  transfer  impedance  (Z^),  as  well  as  pressure,  strain,  and 
wave  propagation  distribution  using  symmetrical  and  asymmetrical  models  of  the  human. 
The  geometric  structures,  such  as  diameters  and  lengths  of  the  airway,  airway  wall  thickness, 
and  branching  patterns,  are  essential  physical  properties  of  the  system  and  thus  significantly 
influence  impedance  and  pressure  distributions.  The  two-port  network  model  based  on  the 
symmetrical  tree  structure  (i.e.,  the  Weibel’s  model  of  the  human  lung  and  the  Yeh’s  model 
of  the  rat  lung)  is  relatively  straightforward  to  realize.  However,  the  resulting  predictions 
may  not  be  entirely  appropriate  since  actual  airways  do  not  branch  symmetrically.  The 
prediction  of  pressure  ratio,  with  forced  oscillation  imposed  at  the  chest  wall  in  the 
symmetrical  airway  network,  did  not  indicate  any  pressure  amplification.  However,  the 
predication  using  the  superposition  model  based  on  the  asymmetrically  branched  structure 
(i.e.,  the  Horsefield’s  model  of  the  human  lung)  showed  that  pressures  in  the  central  airways 
could  be  about  1.6  times  that  of  the  applied  pressure. 

The  computational  difficulty  of  the  asymmetrical  model  was  solved  with  the 
appropriately  designed  algorithms,  the  unique  indexing  representation  of  the  airway 
structure,  division  of  the  airway  network  into  symmetrical  and  asymmetrical  zones,  and 
further  segmentation  of  the  asymmetrical  zone. 

The  Zin  and  Ztr  predictions  of  the  human  lung  models  (both  symmetrical  and 
asymmetrical  airway  networks)  compare  closely  to  actual  measurements  in  adult  humans  for 
frequencies  between  about  2-250  Hz,  and  2-128  Hz,  respectively.  The  predictions  of 
pressure  ratio,  with  forced  oscillation  imposed  at  the  airway  opening,  of  the  symmetrical  and 
asymmetrical  airway  network  compare  qualitatively  to  each  other.  With  forced  oscillation 
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imposed  at  the  chest  wall,  pressure  amplification  (about  1.6  times  that  of  the  input  pressure) 
was  found  at  about  40  Hz  in  the  central  airways  of  the  asymmetrical  airway  network  of  the 
human  lung  model.  However,  this  was  not  found  in  the  symmetrical  human  lung  model.  The 
first  anti-resonant  frequency  was  seen  at  about  90  Hz  in  the  Zin  simulation  of  the  rat  lung 
model,  and  the  first  resonant  frequency  was  seen  at  about  5  kHz.  No  pressure  amplification 
was  found  in  the  rat  lung  model. 

The  reason  for  predicting  the  rat  lung  mechanics  was  to  design  experimental  exposures 
to  study  the  consequences  of  low  frequency  sound  pressure  applied  to  the  body  surface. 
However,  our  results  indicate  that  the  resonances  and  anti-resonances  occur  at  extremely 
high  frequencies  compared  to  the  human  lung. 

We  also  studied  the  consequences  of  airway  constriction,  i.e.,  as  seen  in  asthmatic 
lungs.  The  amplitude  of  resistance  at  the  anti-resonance  decreases  as  the  airways  constrict 
more.  However,  in  the  low  frequency  region,  the  resistance  increases  with  increasing 
constriction.  The  magnitude  of  transfer  impedance  also  increases  as  airways  constrict  more. 

Future  work  in  this  area  of  study  can  include  applying  the  superposition  model  to  the 
Horsefield’s  morphological  data  of  the  dog  lung,  since  it  has  asymmetry  and  in  common 
with  the  Horsefield’s  human  lung  model  as  well  as  the  similar  size. 

The  superposition  model  can  be  also  modified  to  study  the  consequences  of 
heterogeneous  changes  in  airway,  i.e.,  the  constriction  only  occurs  in  left  lobe  of  the  lung 
with  the  right  half  normal. 
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Appendix  A 


1.  Two  network  in  parallel  as  an  analogue  to  two  airway  branches  having 
common  upstream  and  downstream  ends 

When  two  airway  branches  are  connected  in  parallel,  a  common  pressure  node  exists 
at  each  junction,  and  the  continuity  must  be  satisfied  at  each  junction.  From  Fig.  52,  we  can 
see  that: 


Qu  —  Qul  4"  Qu2 

(56) 

Qd  =  Qdl  +  Qd2 

(57) 

Pu  =  Pul  =  Pu2 

(58) 

Pd  =  Pdl  =  Pd2 

(59) 

Figure  52.  Two  networks  are  connected  in  parallel 
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For  network  1, 


(60) 


For  network  2, 


Q2  = 


Qu2 


Qd2 


=  [yj  •  Pi 


(61) 


Compared  to  Eqns  60  and  61,  we  have 


Q  =  Qi +  Q2  =  ( [yi]+[y2]  )  P  =  [y]  P  (62) 


Namely, 


[y]  =  [y']+M 


(63) 


This  can  be  extended  to  any  number  (say,  m)  of  two  ports  connected  in  parallel  as 


(y]=[y']+[y»]+-[y"] 


(64) 
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2.  Transformation  from  Y-matrix  to  T-matrix 


Figure  53.  Two-port  network.  Note  the  signs  of  parameters  depend  on  how  the  flows  and 
pressures  at  two  ends  are  defined. 


The  Y  parameters  for  a  network  (Fig.53)  are 


Thus,  we  have 

Qu  =  yn  •  Pu  +  yi2  •  Pd 
Qd  —  y  21  •  Pu  +  y22  •  Pd 


(65) 


(66) 

(67) 


From  Eq.  67,  P„  can  be  expressed  in  terms  of  Pa  and  Qa  as 

pu  =  -ZE.pd+ J_.Qd  (68) 

y2i  y2i 
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By  substituting  Pu  in  Eq.  66  and  rearranging  it,  we  have 


Pu  =  - 


yn  •  y22  -  yn  •  y2i 
y2i 


Pd  +  yii.Qd  =  _^L.Pd  +  ^.Qd 

y2i  y2i  y2i 


Ay 


yn 


where  Ay  is  the  determinant  of  Y-matrix. 

The  T-parameters  of  the  network  are  defined  as 


(69) 


'Pu' 

"tn 

tl2 

"Pd‘ 

_Q»_ 

t21 

t22 

_Q<«. 

Comparing  Eqns.  68,  69  and  70,  we  have 


tn 

tl2 

y22 

1 

t21 

1 22 

= 

y» 

_  Ay 

y2i 

yn 

y2i 

y2i_ 

(70) 


(71) 


3.  T-matrix  of  tissue  properties 

Consider  the  two-port  network  representation  of  tissue  properties  as  shown  in  Fig. 
54.  Qcw  can  be  expressed  in  terms  of  P  and  Q’  as 


(72) 


100 


PCw  is  obtained  from  Fig.  54 


Figure  54.  Two-port  network  representation  of  tissue  properties 


Thus,  T-matrix  for  tissue  properties  is  obtained  from  Eqns.72  and  73  as 
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4.  Transformation  between  T-matrices 

Consider  Fig.  53  in  Appendix  A.  2.  To  compute  Pd  and  Qd  for  given  P„  and  Qu, 
multiply  Eq.  70  by  the  inverse  of  the  T-matrix,  T'1 , 


T-i  _ 


'til 

tl2 

t22 

—  tl2 

t21 

1 22 

-  t21 

til  _ 

(75) 


Note  that  the  determinant  of  the  T-matrix  is  unity.  We  have 


til  tl2 

-1 

Pu' 

tn  ti2 

-I 

til  tl2 

"Pd' 

1 21  t22 

.Qu. 

r-h 

to 

B 

t21  t22 

_Q<_ 

(76) 


Simplifying  and  reorganizing  Eq.  76  gives  us 


'Pd' 

1 22 

—  tl2~ 

'Pu' 

_Qo_ 

—  t21 

t'l  . 

_Qu. 

(77) 


Therefore,  Pd  and  Qd  can  be  computed  as 

Pd  =  t22  •  Pu  -  tl2  •  Qu  (78) 

Qd  =  — t2i  •  Pu  +  tn  •  Qu  (79) 


102 


REFERENCES: 


1.  Habib  R.  H.,  R.  B.  Chalker,  B.  Suki,  and  A.  C.  Jackson.  “Airway  geometry  and  wall 
mechanical  properties  estimated  from  subglottal  input  impedance  in  humans”,  J.  Appl. 
Physiol.  77(1):  441-451,  1994. 

2.  William  H.  Hayt,  Jr.,  Jack  E.  Kemmerly,  Engineering  Circuit  Analysis.  P76,  1986, 
McGraw-Hill  Book  Company. 

3.  Habib  R.  H.,  B.  Suki,  J.  H.T.  Bates,  and  A.  C.  Jackson.  “Serial  distribution  of  airway 
mechanical  properties  in  dogs:  effects  of  histamine”,  J.  Appl.  Physiol.  77(2):  554-566, 
1994. 

4.  Jackson  A.  C.,  and  K.  R.  Lutchen.  “Physiological  basis  for  resonant  frequencies  in 
respiratory  system  impedances  in  dogs”,  J.  Appl.  Physiol.  70:  1051-1058, 1991. 

5.  DuBois,  A.  B.,  A  W.  Brody,  D.  H.  Lewis,  and  B.  F.  Burgess,  Jr.  “Oscillation  mechanics 
of  the  lungs  and  chest  in  man.”  J.  Appl.  Physiol.  8:  587-594,  1956. 

6.  Lutchen,  K.  R.  and  B.  Suki.  “Understanding  Pulmonary  Mechanics  Using  the  Forced 
Oscillations  Technique.  Emphasis  on  breathing  frequency.”  Bioengineering  Approaches 
to  Pulmonary  Physiology  and  Medicine,  edited  by  Khoo.  Plenum  Press,  New  York, 
1996:  pp.  227-253. 

7.  Lutchen,  K.  R.,  J.  R.  Everett,  A.  C.  Jackson.  “Impact  of  frequency  range  and  input 
impedance  on  airway-tissue  separation  implied  from  transfer  impedance.”  J.  Appl. 
Physiol.  74(3):  1089-1099,  1993. 


103 


8.  Peslin,  .R.  and  JJ  Fredberg.  “Oscillation  Mechanics  of  the  respiratory  system.” 
Handbook  of  Physiology.  The  Respiratory  System.  Mechanics  of  Breathing.  Bethesda, 
MD:  Am.  Physiol.  Soc.,  1986,  Section  3,  Volume  III,  Part  1,  Chapter  11,  pp.  145-176. 

9.  Sullivan,  A.  “Low  and  high  frequency  transfer  impedance  measurements  in  humans.” 
Master ’s  Thesis.  Boston  University,  Department  of  Biomedical  Engineering,  Boston, 
MA,  1995. 

10.  Helal  M.  A.,  K.  C.  Watts,  and  S.N.  Sarwal.  “Theoretical  model  for  assessing 
haemodynamics  in  arterial  networks  which  include  bypass  grafts.”  Med.  &  Biol.  Eng.  & 
Comput.,  1990,  28,  pp.Yen.  465-473 

11.  Weibel,  E.  R.  “Morphometries  of  the  lung.”  Handbook  of  Physiology.  Respiratory,  sect. 
3,  vol.  I.,  edited  by  W.  O.  Fenn  and  H.  Rahn.  Washington,  DC:  Am.  Physiol.  Soc.,  1964, 
pp.  299-305 

12.  Horsefield,  K.,  and  G.  Cumming.  “Angles  of  branching  and  diameters  of  branches  in  the 
human  bronchial  tree.”  Bull.  Math.  Biophys.  29:  245-259,  1967 

13.  Personal  communication  with  Prof.  Suki,  BME  department,  Boston  University 

14.  Nedwell,  J.  R.,  and  S.  J.  Parvin.  “The  effects  of  low  frequency  sonar  transmissions  on 
divers  and  ichthyofauna:  literature  survey  and  initial  experimental  results.”  DRA(AWL) 
TM93721,  Sept.  1993 

15.  Suki,  B.,  R.  H.  Habib,  and  A.  C.  Jackson.  “Wave  propagation,  input  impedance,  and  wall 
mechanics  of  the  calf  trachea  from  16  to  1600  Ha.”  J.  Appl.  Physiol.  75:  2755-2766, 
1993 


104 


16.  Peslin,  R.  and  J.  J.  Fredberg.  “Morphometry  of  airways”  The  Respiratory  System. 
Mechanics  of  Breathing.  Bethesda,  MD:  Am.  Physiol.  Soc.,  1986,  Section  3,  Volume 
m,  Part  1,  Chapter  7,  pp.  75-88. 

17.  Hayt,  W.  H.  and  J.  E.  Kemmerly.  “Two-port  networks.”  Engineering  Circuit  Analysis, 
McGraw-Hill  Book  Company,  1986,  Chapter  16,  pp.  458-485 

18.  Comroe,  J.  H..  “Physiological  diagnosis.”  Physiology  of  Respiration,  Year  Book  Medical 
Publishers,  Inc.,  1966,  Chapter  21,  pp  218 

19.  Lambert,  R.  K.,  A.  R.  Frey,  A.  B.  Coppens,  and  J.  V.  Sanders.  “A  method  for  estimating 
the  Young’s  modulus  of  complete  tracheal  cartilage  rings.”  J.  Appl.  Physiol.  70:  1152- 
1159,  1991 

20.  Jackson,  A.  C.,  M.  Tabrizi,  M.  I.  KotlikofF,  and  J.  R.  Voss.  “Airway  pressures  in  an 
asymmetrical  branched  airway  model  of  the  dog  respiratory  system.”  J.  Appl.  Physiol .: 
Respirat.  Environ.  Exercise  Physiol.  57(4):  1222-1230,  1984 

21.  Guelke,  R.  W.,  and  A.  E.  Bunn.  “Transmission  line  theory  applied  to  sound  wave 
propagation  in  tubes  with  compliant  walls.”  Acoustica  48:  101-106,  1981 

22.  Dorkin,  H.  L.,  K.  R.  Lutchen,  and  A.  C.  Jackson.  “Human  respiratory  input  impedance 
from  4  to  200  Hz:  physiological  and  modeling  considerations.”  J.  Appl.  Physiol.  64(2): 
823-831,  1988 

23.  Jackson,  A.  C.,  and  J.  W.  Watson.  “Oscillatory  mechanics  of  the  respiratory  system  in 
normal  rats.”  Respiration  Physiol.  48:  309-322, 1982 

24.  Takezawa  J.,  F.  J.,  Miller,  and  J.  O’Neil.  “Single-breath  diffusing  capacity  and  lung 
volumes  in  small  laboratory  mammals.”  J.  Appl.  Physiol .:  Respirat.  Environ.  Exercise 
Physiol.  48(6):  1052-1059,  1980 


105 


4 


25.  Zwart,  A.  and  E.  M.  Hessel.  “Oscillatory  mechanics  of  the  mouse  respiratory  system.’ 
Eur.  Respir.  Rev.  4:  19,  182-186,  1994 

26.  Habib,  R.  H.  “Airway  geometry  and  wall  mechanical  properties  from  respiratory  acoustic 
impedance.”  Ph.D.'s  dissertation.  Boston  University,  Department  of  Biomedical 
Engineering,  Boston,  MA,  1989. 


106 


